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The ground-state phase diagram of a frustrated S = 1 Heisenberg spin chain with uniaxial 
single-ion-type anisotropy is studied using the infinite-system density-matrix renormalization- 
group method. Particular attention is paid to "chiral" phases, in which the long-range order 
parameter is a vector chirality. It is found that, in addition to the Haldane, large-D (LD), 
and double-Haldane phases, the phase diagram includes three different chiral phases, i.e., the 
gapless chiral, chiral-Haldane, and chiral-LD phases. These chiral phases are distinguished from 
each other by the energy gap and the string long-range order. 
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ordering, density-matrix renormalization-group method 

§1. Introduction 

The effect of geometrical frustration in quantum spin systems has attracted considerable attention 
for many years. Frustration generally suppresses the tendency towards the classical Neel ordering 
and leads to a wide variety of exotic phenomena such as a spontaneous breaking of the translational 
symmetry in a one-dimensional (ID) system and spin-liquid states in 2D and 3D systems. Among 
them, the possibility of a novel "chiral" phase in ID frustrated quantum spin systems has been 
studied intensively in recent years. This chiral phase is characterized by the absence 

of the helical long-range order (LRO), 

miq) = ^Y.Sl^"''=^ (forV<7), (1.1) 

and the non-zero value of the z-component of the total vector chirality,& 

O« = ;^E^'^0' (1-2) 
Ki = SfSl, - SfSf_,, = [Si X Si+il , 
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where is a spin-5 operator at site / and L is the system size assumed to be even throughout this 
paper. The chiral phase breaks only the parity symmetry spontaneously with preserving the time- 
reversal and translational symmetries. Using the bosonization method, Nersesyan et al. predicted 
the appearance of the chiral phase in the frustrated antiferromagnetic 5 = 1/2 chain with nearest- 
neighbor and next-nearest-neighbor couplings of easy-plane anisotropy.0^ From the analytica]@'& 
and numericafli'S studies, it has been shown that the chiral phase appears in the frustrated chain 
with the anisotropic couplings in general spin-S" cases. The ground-state phase diagrams for the 
S = 1/2, 1, 3/2, and 2 cases have been determined numerically.!'!^ 

Another interesting observation in the studiesi'SHi-* is that there exist two different chiral phases 
in the integer-S cases. One of them is the "chiral-Haldane" phase where the chiral LRO coexists 
with the string LRo!'0) 

Ostr = ^T.^^V St, (1.3) 



which characterizes the well-known Haldane phasejHlllHP and the spin correlation function decays 
exponentially with a finite energy gap. The other is the "gapless chiral" phase where the chiral LRO 
exists and the spin and string correlation functions decay algebraically with gapless excitations. 
It has been shown that for S = 1 and 2 the gapless chiral phase exists in a broad region of the 
ground-state phase diagrams while the chiral-Haldane phase exists in a narrow region between the 
Haldane and gapless chiral phases. (See Fig. 1 in ref. ^ and Fig. 9 in ref. ^.) For S = 1/2 and 3/2, 
on the other hand, the gapless chiral phase exists as in the integer-5 cases while the chiral phase 
with a finite energy gap has not been identified. 

As the mechanism generating the chiral phases, easy-plane anisotropy plays an essential role 
as well as frustration does. The anisotropy lowers the symmetry of the system from SU{2) to 
U{1) X Z2 and induces a two- fold discrete degeneracy according as the right- and left-handed 
chirality. In fact, the chiral phases do not appear in the isotropic case where the discrete chiral 
degeneracy no longer exists. From the experimental viewpoint, the anisotropy in real materials 
commonly takes the form of the uniaxial single-ion-type. It is therefore important to study the 
effect of the single-ion-type anisotropy on the chiral phases, for the sake of exploring experimental 
implications of the theoretical results and stimulating an investigation into the chiral phases in real 
materials. Nevertheless, most of the theoretical findings on the chiral phases mentioned above have 
been obtained for the case of the anisotropic exchange couplings. Kolezhuki^ studied recently the 
frustrated integer-5 spin chain with the single-ion-type anisotropy by means of a large-5 approach 
and presented a schematic phase diagram including the gapless chiral phase and a chiral phase with 
a finite energy gap. However, systematic studies on the chiral phases in the case of single-ion-type 
anisotropy have been scarce so far. 

The aim of this paper is to investigate the ground-state properties of the frustrated 5 = 1 
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Heisenberg chain with the uniaxial single-ion-type anisotropy. The model Hamiltonian is given by 

^= E \jpj:Si-SiA+Dj2{Sff. (1.4) 

p=l,2 [ I ) I 

I concentrate on the case of the antiferromagnetic couplings ( Ji > and J2 > 0) and the easy-plane 
anisotropy {D > 0). Hereafter I denote j = J2/J1 and d = D/Ji. The ground-state properties 
of the model ( |1.4| ) have been studied for some limiting cases. For j = the phase transition 
between the Haldane and large-D (LD) phases have been studied intensively. •l3 It has been 

found that as d increases the system undergoes a continuous phase transition from the Haldane 
phase to the LD phase at d = 1.000 it O.OOlIiI) Meanwhile, it has been shown for d = that as 
j increases the system undergoes a first order transition at j ~ 0.747 from the Haldane phase to 
the "double Haldane" (DH) phase, I Bill) in which the next-nearest-neighbor coupling J2 becomes 
dominant so that the system is described as two Haldane subchains weakly coupled by the inter- 
subchain coupling Ji. The string LRO vanishes discontinuously at the transition. Furthermore, the 
appearance of the gapless and gapped chiral phases has been predicted by the large-S" approach.!^ 



In the present work, I determine numerically the ground-state phase diagram of the model (1^) by 
employing the same method as the one used in our previous worksiBSi^ I calculate the appropriate 
correlation functions associated with the order parameters characterizing each phase using the 
density-matrix renormalization group (DMRG) methodjllB and then, analyze their long-distance 
behaviors. From the obtained phase diagram, we find that the gapless chiral phase appears in a 
broad region of the j-d plane while the chiral-Haldane phase appears in a narrow region between 
the Haldane and gapless chiral phases, which are essentially the same as those in the case of the 
anisotropic exchange couplings. It is found further that, in addition to the gapless chiral and chiral- 
Haldane phases, another novel chiral phase exists between the gapless chiral and LD phases. This 
"chiral-LD" phase is characterized by the chiral LRO and the exponential decay of both the spin 
and string correlation functions. 

The paper is organized as follows. In §2, I introduce the chiral, string, and spin correlation 
functions associated with each order parameter and explain the numerical method. The obtained 
phase diagram and numerical data are presented in §3. The results are summarized in §4. 

§2. Correlation Functions and Numerical Method 

The method used to determine the phase diagram is essentially the same as those used in our 
previous works.Bi'S I have calculated the two-point chiral, string, and spin correlation functions 
defined by 

c^ii,l') = {Km,), (2.1) 

Cstr(/,/') = (5f exp lmJ2sn SP), (2.2) 
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C^,{l,l') = {SfSf,), (2.3) 

which are associated with the order parameters (1.2), (|1.3|), and ([L.lj), respectively. The notation 
(• • •) represents the expectation value in the lowest energy state in the subspace of S^^^^^ — J2i Sf = 
with even parity. It has been checked by the exact-diagonalization calculation that the ground state 
of the chain indeed belongs to the subspace. The calculation has been performed for various fixed 
values of d (j) with varying j (d). Then, I have estimated the transition points jc (dc) by examining 
the dependence of the correlation functions on r = |/ — at long distance. 

I have employed the infinite-system DMRG algorithm originally introduced by Whitellli^) and 
accelerated it by making use of the recursion relation proposed by Nishino and Okunishi.0'ii) The 
number of kept states m is up to 350 if not otherwise mentioned. The convergence of the data 
with respect to m has been checked by increasing m consecutively. Since the truncation error of 
the DMRG calculation increases dramatically as j increases, the calculation is limited to rather 
small j. However, I consider that the region of j treated in this work, j < 1, is enough to study 
the properties of the chiral phases which we wish to know. In the calculations, the open boundary 
condition has been imposed for the sake of the best performance of the algorithm. Accordingly, 
in order to avoid the unwanted boundary effect, the two-point correlation functions have been 
calculated for the sites / and /' near the center of the chain, i.e., / = Iq — r/2 and I' = Iq + r/2 
where Iq = L/2 for even r and Iq = {L + 1) /2 for odd r. By checking the convergence of the data 
with increasing L typically up to 2000, I have confirmed that the data are free from the effect of 
the open boundaries. 

§3. Phase Diagram 

By examining the long-distance behaviors of the chiral, string, and spin correlation functions 
introduced in the previous section, I have determined the ground-state phase diagram of the model 
( |1.4| ) in the j-d plane. The obtained phase diagram shown in Fig. 1 includes six different phases, 
i.e., the Haldane, LD, DH, gapless chiral, chiral-Haldane, and chiral-LD phases. The long-distance 
behaviors of the correlation functions in each of the phases are summarized in Table I. We can 
see in Fig.l that the gapless chiral phase appears in a broad region of the phase diagram while 
the chiral-Haldane phase appears between the Haldane and gapless chiral phases. Further we 
find the chiral-LD phase, where the chiral LRO exists while the spin and string correlation decay 
exponentially, between the LD and gapless chiral phases. Two multi-critical points, at least, are 
observed: One of them is at {jMi-,dMi) — (0.51,0.76) among the Haldane, LD, gapless chiral, 
chiral-Haldane, and chiral-LD phases while the other is at (jm2)'^M2) — (0.75,0.15) among the 
Haldane, gapless chiral, chiral-Haldane and DH phases. The numerical data for each transition line 
are discussed in the following subsections. 
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3.1 Haldane-chiral transition 

First let us consider the phase transitions between the Haldane and chiral phases. Figures 2(a)- 
(c) show the data of the chiral, string, and spin correlation functions in log-log plots for d = 0.5 and 
several typical values of j. The data of the spin correlation shown in Fig. 2 (c) are divided by the 
leading oscillating factor cos{Qr) where Q is the wavenumber characterizing the incommensurate 
oscillation of (r) in real space. It can be clearly seen in Fig. 2 (a) that the data of the chiral 
correlation function Ci^{r) are bent upward for j > jd — 0.635 suggesting a finite chiral LRO while 
they are bent downward for j < jci suggesting the absence of the chiral ordering. The transition 
point where the chiral LRO sets in is estimated to be jd = 0.635ib0.005. As shown in Fig. 2 (b), the 
string correlation function exhibits a finite LRO for j < jc2 — 0.670 whereas it decays algebraically 
for j > ^£2-0^ Meanwhile, the data of the spin correlation function shown in Fig. 2 (c) are bent 
downward for j < jc2 suggesting an exponential decay with a finite energy gap while they exhibit 
a linear behavior for j > jc2 suggesting a power-law decay with gapless excitations. The transition 
point jc2 where the string LRO and the energy gap vanish is estimated to be jc2 = 0.670 =b 0.010. 
It should be noticed that the estimate of jc2 is distinctly larger than that of jd- Accordingly, 
the chiral-Haldane phase where the chiral and string LROs coexist and the spin correlation decays 
exponentially exists in a narrow but finite region between the Haldane and gapless chiral phases. 
The existence of the chiral-Haldane phase can be clearly seen from the behavior of the correlation 
functions at j = 0.640 and 0.660 which lie between jd and jc2- Thus, it is concluded that the system 
undergoes two successive transitions as j increases, first at j = jd from the Haldane phase to the 
chiral-Haldane phase, and then, at j = jc2 to the gapless chiral phase. The results are essentially 
the same as the ones observed in the frustrated chain with anisotropic exchange couplings .Bi^ 

Performing calculations in the same manner for various fixed d (j) with varying j (d), I have 
determined the transition points jd {dd) and jc2 {dc2) plotted in Fig. 1. The figure may sug- 
gest that the transition lines connecting the estimated points run from the multi-critical point 
{jMi,dMi) to the other multi-critical point (jA/2;C?M2)- The chiral-Haldane phase is identified for 
d = 0.7,0.6,0.5,0.4, and 0.3 while it is not for d = 0.2 within the calculation with varying j at in- 
tervals of 0.002. However, I note that it is difficult to exclude the possibility that the chiral-Haldane 
phase exists in a too narrow region to be detected. Thus it is not clear whether the transition lines 
jd and jc2 merge at {jM2-,dM2) or at another multi-critical point '^M2) with d'j^^g > 0.2. In the 
latter case, a transition line j'd between the Haldane and gapless chiral phases connects {jM2^d-M2) 
and {jM2,dM2)- For d = 0.1 and 0.05, on the other hand, the chiral LRO is not observed. In 
these cases, it is also found that the string LRO vanishes discontinuously at the transition point 
and the spin correlation function exhibits an exponential decay for all calculated values of j. The j 
dependence of the string correlation function at r — > oo for d = 0.1 is shown in Fig. 3 together with 
the spin-correlation length ^ estimated by fitting the data of (r) to the form of an exponential 
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decay, Cf(r) = j4 cos((5r) exp(— r/^), where j4 is a numerical constant. These behaviors of the 
string and spin correlation functions for d = 0.1 and 0.05 suggest a first-order transition between 
two gapful phases, i.e., the Haldane and DH phases. Prom the observations, it is concluded that 
between d = 0.2 and d = 0.1 there is a multi-critical point (jm2, c^M2) where the transition lines jd 
and jc2 or the transition line j'^i merges with the DH-chiral transition line discussed in §3.4. The 
first-order transition line between the Haldane and DH phases smoothly connects the multi-critical 
point (jm2) dM2) to the transition point in the isotropic case {d = 0), jc2 — 0.747, estimated in the 
previous works. 

3.2 LD-chiral transition 

We next consider the phase transitions between the LD and chiral phases. The data of the chiral, 
string, and spin correlation functions for j = 0.8 and several typical values of d are shown in Figs. 
4 (a)-(c) in log-log plots. As shown in Fig. 4 (a), the chiral correlation function Cfi{r) exhibits a 
finite LRO for d < dc3 ~ 1.09 while it decays exponentially for d > dc^. The transition point dcs 
where the chiral LRO vanishes is thereby estimated as dcs = 1-09 it 0.01. Meanwhile, the string 
and spin correlations decay exponentially for d > dd — 0.80 suggesting a finite energy gap whereas 
they decay algebraically for d < dcA suggesting gapless excitations.© The transition point dcA is 
estimated to be d^ = 0.80 ± 0.10. Unfortunately, the numerical error of the estimate of dcA is 
quite large due to the truncation error of the DMRG calculation, which tends to underestimate the 
string and spin correlation functions. Nevertheless, from the behaviors of the correlation functions 
for 0.94 < d < dcz, it can be clearly seen that there exists a intermediate region with a finite chiral 
LRO and a finite energy gap. In order to confirm the appearance of the intermediate phase, I have 
performed a high-precision calculation with m up to 500 for d = 1.00 and 0.94. The data shown 
in Fig. 4, which are almost free from the truncation error, clearly suggest that the parameter 
points d = 1.00 and 0.94 indeed belong to the chiral phase with a finite energy gap. It is therefore 
concluded that the estimate of dcA is distinctly smaller than that of dcs and the system undergoes 
two successive transitions as d decreases, first sX d = d^z from the LD phase to the "chiral-LD" phase 
characterized by a finite chiral LRO and an exponential decay of the string and spin correlation 
functions, and then at d = dd to the gapless chiral phase. Here it should be noticed that the 
chiral-LD phase is distinct from the chiral-Haldane phase with respect to the absence of the string 
LRO. Hence the chiral-LD phase can be regarded as a new type of the gapped chiral phase. 

By the same calculation for various fixed j (d), the transition points (ics (jcs) and dd {jci) are 
estimated. In particular, the appearance of the chiral-LD phase has been confirmed by the high- 
precision calculations with m up to 500 for the points plotted in Fig. 1 by crosses. As seen in Fig. 
1, the chiral-LD phase appears in a finite region between the LD and gapless chiral phases. The 
transition lines (ic3 and dcA merge with the Haldane-chiral transition lines jd and jc2 around the 
multi-critical point (jmi, dMi) into the Haldane-LD transition line discussed in the next subsection. 
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3.3 Haldane-LD transition 

In this subsection, we see the results for the phase transition between the Haldane and LD phases. 
The phase transition has been studied in detail in the case of no frustration (j = 0). i'BBIll) It 
has been pointed out that there occurs a continuous phase transition accompanying the vanishing 
of the energy gap at d = 1.000 ± 0.001. In Figs. 5 (a) and (b), I show the data of the string and 
spin correlation functions for j = 0.2 and several typical values of d. Here the m convergence of 
the data has almost been achieved even at m = 180. We can see in Fig. 5 (a) that the string 
correlation function exhibits a finite LRO for d < dc^ — 0.90 while it exhibits an exponential decay 
for d > dc5. In the meantime, the spin correlation function exhibits an exponential decay for 
d < 0.89 and d > 0.91 in accordance with the existence of a finite energy gap in both phases while 
it decays algebraically at d = 0.90 suggesting the vanishing of the energy gap at the transition. 
The transition point between the Haldane and LD phases for j = 0.2 is thereby estimated to be 
dc5 = 0.90 lb 0.01. The transition points dc5 estimated in the same way are plotted in Fig. 1. The 
transition line connecting the estimates runs from the point = 0.96 it 0.02 at j = 0, which is 
slightly smaller than the previous estimate dc^iJ = 0) = 1.000 it 0.001, to the multi-critical point 

{jMl,dMi). 

3.4 DH-chiral transition 

Here we consider the phase transition between the DH and chiral phases. Figure 6 (a) shows the 
chiral correlation function (^^(r) for d = 0.3 and several typical values of j around the DH-chiral 
transition. As shown in the figure, Ck{t) exhibits a finite LRO for j < jcQ — 0.820 while it decays 
exponentially for j > jcQ. Thus, the transition point where the chiral LRO vanishes is estimated to 
be jcQ = 0.820 lb 0.010. Meanwhile, the data of the spin correlation function for d = 0.3 are shown 
in Fig. 6 (b) where the data are divided by the oscillating factor cos{Qr). It can be seen in the 
figure that as j increases the spin correlation changes its behavior from an algebraic decay to an 
exponential decay at j = jc7 — 0.780. Although the truncation error of the data prevents us from 
carrying out the precise estimation, the transition point is estimated to be jc7 = 0.780 it 0.040. 

Here arises a question: Whether jc7 is equal to or smaller than jcs ? If jcv is equal to jc6, there 
occurs only one phase transition at j = Jcq = jc7 between the gapless chiral phase and the DH phase. 
On the other hand, if jc7 is smaller than jcQ, the system undergoes two successive transitions as j 
increases, first at j = jcQ from the gapless chiral phase to an intermediate "chiral-DH" phase with a 
finite chiral LRO and a finite energy gap, and then at j = jc7 to the DH phase. Unfortunately, the 
estimates of the transition points, particularly that of jc7, are not accurate enough to determine 
which of the cases is realized. Here I note that the data of the spin correlation for j = 0.800 shown 
in Fig. 6 (b) are calculated with m up to 500 and still suffer from the truncation error which is not 
negligible. This means that it is unpromising to achieve the sufficiently accurate estimate of jc7 in 
the manner used in this work. The difficulty in identifying the intermediate phase in this case can 
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be ascribed partly to the fact that no order parameter characterizing the DH phase has been found 
so farB and partly to the fact that the intermediate phase, if any, appears only in a narrow region 
which can be concealed behind the error bar of jcj. Hence, the situation here is different from that 
of the Haldane-chiral transition where we can utilize the string order parameter to estimate jc2 and 
that of the LD-chiral transition where the region of the chiral-LD phase is so broad that we can 
identify the chiral-LD phase even with rather rough estimates of dcA- Since it is a subtle problem 
to estimate the point where the behavior of the spin correlation changes, it might be needed for 
the settlement of the question to introduce an order parameter characterizing the DH phase. This 
problem remains open for future studies. 

In the same manner, I have estimated the transition points jcG and jc7 for d = 0.2 and d = 0.4. 
The chiral-DH phase has not been identified also in these cases. As shown in Fig. 1, the transition 
lines connecting the estimates seem to emerge from the multi-critical point (ij\/2, dM2)- The values 
of jcQ and jc7 become larger as d becomes larger suggesting that the DH phase becomes stable as 
j increases. 

3.5 Multi- critical point 

Lastly, I touch upon the region of the phase diagram around the multi-critical point (Jmi, dMi)- 
It can be seen in the obtained phase diagram Fig. 1 that the five transition lines run into a narrow 
region around {Jmi, dMi)- It might be useful to classify them into three groups, i.e., the transition 
lines jci and dcs where the chiral LRO sets in, the lines jc2 and (ic4 between gapful and gapless 
phases, and the line dc^ between the Haldane and LD phases. Here I note that the mechanism 
generating the chiral LRO is essentially distinct from the one generating the energy gapHi^ and, 
to my knowledge, no definite relation between them has been known. Thus, the behavior of the 
chiral transition lines {jd and dcs) and the "gapful-gapless" lines (jc2 and dc4,) around the multi- 
critical point is still nontrivial. Looking at Fig. 1, two possibilities seem to be left: (i) The five 
transition lines merge at one multi-critical point (jAfi,(iMi)- (h) The transition lines jd and dcs 
merge with the Haldane-LD transition line dc5 at {Jmi, dMi) while the lines jc2 and (ic4 merge with 
the transition line d'^^ between the chiral-Haldane and chiral-LD phases at another multi-critical 
point Oa/I'^^mi)- These scenarios are illustrated in Fig. 7 (a) and (b), respectively. Because of the 
rather large error of the estimates of (ic4 mentioned in §3.2, we can not decide which of the cases is 
realized from the present numerical results. More detailed studies might be needed to identify the 
universality class of the multi-critical point(s) and transition lines. 

§4. Summary 

In this paper, the ground-state properties of the frustrated S = 1 Heisenberg spin chain with 
uniaxial single-ion- type anisotropy (eq. ( |1.4D ) have been investigated numerically by the infinite- 
system DMRG method. By calculating the chiral, string, and spin correlation functions and analyz- 
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ing their long-distance behaviors, I have determined the ground-state phase diagram for < j < 1 
and d > (Fig. 1). It has turned out that there exist six different phases, namely, the Haldane, 
LD, DH, gapless chiral, chiral-Haldane, and chiral-LD phases. The gapless chiral phase appears 
in a broad region of the phase diagram while the chiral-Haldane phase appears in the narrow but 
finite region between the Haldane and gapless chiral phases; these results are essentially the same 
as those found in the frustrated S = 1 chain with anisotropic exchange couplings.@i) Meanwhile, 
the chiral-LD phase has been found between the LD and gapless chiral phases. Since the chiral-LD 
phase is characterized by a finite chiral LRO, the absence of the string LRO, and the short-range 
spin correlation function, it is distinct from either the gapless chiral or the chiral-Haldane phase 
and can be regarded as a new type of the gapped chiral phases. The question about the existence 
of the possible chiral-DH phase remains open. 

After the theoretical finding of the chiral phases, one of the most challenging tasks is the ex- 
perimental observation of the phases in real materials. As seen in the obtained phase diagram, 
the anisotropy (i>0.15 and the next-nearest neighbor coupling j within a suitable range around 
j ~ 0.75 are required for a material realizing the chiral phases. These values of d and j are re- 
alistic so that there must be a good chance to find out a material in the parameter region. If 
such a material is prepared, the measurement of the vector chirality (1.2) is, in principle, possible 
by using polarized neutrons. El' 00 I hope that the results of the present paper stimulate further 
experimental studies. 
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Fig. 1. The ground-state phase diagram of the model ([l.4[). The circles, squares, and diamonds represent the 
points where the chiral LRO sets in, the points where the string LRO vanishes, and the points where the behavior 
of the spin and string correlations changes from an exponential to an algebraic decay. The multi-critical points 
(jMiydMi) and (jAi2,dA/2) are indicated as "Ml" and "M2", respectively. The crosses represent the points for 
which the calculations with m — 500 are performed. The lines are to guide the eye. 
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Fig. 2. The r dependence of the correlation functions for d = 0.5: (a) chiral correlation (^^(r); (b) string correlation 
— C'Btr(r); (c) spin correlation Cg{r) divided by the oscillating factor cos(<5r). The number of kept states is m = 350. 
To illustrate the m dependence, the data with m = 260 and 300 are indicated by crosses: for j = 0.620 and 0.630 

in (a); for j — 0.680 and 0.700 in (b). In other cases, the truncation errors arc smaller than the symbols. The data 
of CK{r) for j = 0.700 in (a), which are almost the same as those for j = 0.680, are omitted for clarity. 



Fig. 3. The j dependence of the string correlation function —Csti{r) at r ^ oo (solid circle) and the inverse spin- 
correlation length (open circle) for d = 0.1. The dotted line represents the transition point. 



Fig. 4. The r dependence of the correlation functions for j = 0.8: (a) chiral correlation CK(r); (b) string correlation 
— C'str(r); (c) spin correlation (r) divided by the oscillating factor cos(Qr). The number of kept states is m = 500 
for d = 0.94 and 1.00 while it is m = 350 for the other cases. To illustrate the m dependence, the data with m = 260 
and 300 for d = 0.80, and 0.70 are indicated in (b) by crosses. In other cases, the truncation errors are smaller 
than the symbols. 



Fig. 5. The r dependence of the correlation functions for j — 0.2: (a) string correlation —Catrir); (b) spin correlation 
(— l)''Cs (r). The number of kept states is m = 180. The truncation errors are smaller than the symbols. 



Fig. 6. The r dependence of the correlation functions for d — 0.3: (a) chiral correlation CK{r); (b) spin correlation 
Cs (r) divided by the oscillating factor cos{Qr). The number of kept states is m = 500 for j = 0.800 and m = 350 
for other values of j. To illustrate the m dependence, the data with smaller m are indicated by crosses: the data 
with m = 260 and 300 for j = 0.820 in (a); the data with m = 350 and 400 for j = 0.800 and with m = 260 and 
300 for j = 0.760 in (b). In other cases, the truncation errors are smaller than the symbols. 



Fig. 7. The possible schematic pictures of the phase diagram around the multi-critical point {jMi,dMi)- Figures 
(a) and (b) correspond to the cases (i) and (ii) in text, respectively. 
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Table I. The long-distance behaviors of the spin, chiral, and string correlation functions in each of the phases. The 
words "expo." and "power" mean an exponential decay and a power- law decay, respectively. 
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expo. 



12 




Fig.1 , T.Hikihara 



(a) chiral correlation 



d = 0.5 




(b) string correlation 



d = 0.5 



10" 



O 



* * * ••••—mmmm 

° O OOOOi 



10" 



= 0.620 
= 0.630 
= 0.640 
= 0.660 
= 0.680 
= 0.700 

10 




(c) spin correlation 



d = 0.5 




= 0.620 
= 0.630^ 
= 0.640 
= 0.660^ 
= 0.680" 
= 0.700 



■Jc1 
-jc2 



Fig.2, T.Hikihara 




0.05 



Fig. 3, T.Hikihara 



10' 



(a) chiral correlation 



O 



10" 



j = 0.8 



m :d = 1.10 




o : d = 1 .08 
■ :d = 1.06 
□ :d = 1.00 
A : d = 0.94 
A : d = 0.80 ^ clc4 



(b) String correlation 



10-f 'B^iiA 



(3 



10" 



10" 



• :d = 1.10^H 
o : d = 1 .00 
■ : d = 0.94 
□ : d = 0.80 ^ dc4 
A : d = 0.70 



j = 0.8 




10 



10"= 



spin correlation 




Fig.4, T.Hikihara 



(a) string correlation j = 0.2 





Fig. 5, T.Hikihara 



(a) chiral correlation 



d = 0.3 





Fig. 6, T.Hikihara 




chiral >— ^ 
Haldane (Ml) 



w gapless 
\ chiral 



Fig. 7, T.Hikihara 



